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ABSTRACT 

Dwarf stars are believed to have small protostar disk where planets may grow up. During the 
planet formation stage, embryos undergoing type I migration are expected to be stalled at inner edge 
of magnetic inactive disk (ocrit ~ 0.2 — 0.3 AU). This mechanism makes the location around Ocrit a 
sweet spot of forming planets. Especially, Ocrit of dwarf stars with masses ~ 0.5 Mq is roughly inside 
the habitable zone of the system. In this paper we study the formation of habitable planets due to this 
mechanism with a model system OGLE-06-109L. It has a 0.51 Mq dwarf star with two giant planets 
in 2.3 and 4.6 AU observed by microlensing. We model the embryos undergoing type I migration 
in the gas disk with a constant disk accretion rate (M). Giant planets in outside orbits affect the 
formation of habitable planets through secular perturbations at the early stage and secular resonance 
at the later stage. We find that the existence and the masses of the habitable planets in OGLE-06-10L 
system depend on both M and the speed of type I migration. If planets formed earlier so that M is 
larger ('~ 10~^Mq yr~^), terrestrial planet can not be survived unless the type I migration rate is an 
order of magnitude less. If planets formed later so that M is smaller (~ 10~^Mo yr~^), single and 
high mass terrestrial planets with high water contents (~ 5%) will be formed by inward migration of 
outer planet cores. A slower speed migration will result in several planets by collisions of embryos, 
thus their water contents are low (~ 2%). Mean motion resonances or apsidal resonances among 
planets may be observed if multiple planets survived in the inner system. 

Subject headings: (stars:) planetary systems: formation - solar system: formation - stars: individual 
(OGLE-06-109L) 
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1. INTRODUCTION 

One of the interesting problems for planetary dynam- 
ics is whether there exist stable orbits in the habitable 
zone (HZ) of main sequence stars. Earth-size planets in 
HZ possibly have similar compositions especially liquid 
water, which is important for life survival (Kasting et 
al. 1993; Hinse et al. 2008). Currently, radial veloc- 
ity technique is not sensitive enough to detect terrestrial 
planets in HZ of ~ 1 AU around solar type stars. How- 
ever, dwarf stars with masses ~ 0.5 Mq have very close 
habitable zones 0.2 — 0.3 AU) due to their low stellar 
luminosity, thus they are good candidates for terrestrial 
planets detection in HZ. 

According to the simulations of terrestrial planet for- 
mation in exoplanet systems, in site formation of Earth- 
size planets in close-in orbits requires a very massive cir- 
cumstellar disk, e.g., ~ 30 times of minimum solar neb- 
ular (Raymond et al. 2008) to form a ^ 10 M^ planet 
around GJ 581. Thus formation of Earth-like planets 
in close-in orbits requires either keeping adding building 
blocks into inner orbits, or massive cores formed in out- 
side orbits and migrating inward under, e.g., the type I 
migration. 

If giant planets already formed in outside orbits, per- 
turbations from giant planets that undergo type II migra- 
tions may be helpful to the final assembly of terrestrial 
planets in close-in orbits. During the inward migration 
of giant planets embedded in the viscous gas disk, the 
locations of inner mean motion resonances (MMR) with 
the gas giants , mainly 2:1 MMR, will sweep through 
the inner disk, trap and shepherd the embryos, and ex- 



cite their eccentricities, which results in the mergers of 
isolated embryos and the formation of close-in planets 
(Zhou et al. 2005; Fogg & Nelson 2005; Raymond et al. 
2006). During the stage of gas disk depletion, secular res- 
onances caused by disk potential and outer giant planets 
will sweep into the inner region, which may excite the ec- 
centricities of the embryos and results in the merger and 
formation of exoplanets (Nagasawa et al. 2003; Zhou et 
al. 2005). 

In the above mentioned scenarios, the type I migra- 
tion of embryos is not considered. Due to the fast rate 
of type I migration, the sweeping of 2:1 MMR (caused 
by type II migration of outside giant planets) and secu- 
lar resonances may be lagged the fast inward migration 
of embryos, thus their effects will be greatly reduced. In 
that case, the embryos will migrate unless there are some 
mechanisms to shop it. Although in Fogg & Nelson 2007 
they considered the type I migration in hot Jupiter sys- 
tems, the formation of terrestrial planets in close-in or- 
bits under the perturbations of giant planets in the outer 
region need to be re-investigated to see whether secular 
perturbations and secular resonances are effective. 

Recent hydrodynamical simulations indicate that a lo- 
cation of density maximu m in the gas disk ca n stop the 
type I inward migration ('Mass et et al. 20061 ). During 
the evolution of the protoplanct, the exchange of angu- 
lar momentum between the protoplanet and the gas disk 
generates a cor otation torque on the protoplanet wit h 
the expression (IWard 19911 : iMasset 200l1: iMasset 200"1) . 
Fc oc 'E,gd\n{Tig/ B)/dhia, where Eg is the surface den- 
sity of the gas disk, B is the Oort constant and B oc a^"^/^ 
in a near Keplerian motion. If the surface density can 
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jump over 3 ~ 5 gas disk thickness, the positive corota- 
tion torque will overcome the negative Lindblad torque 
which casues the inward type I migration, making the 
location of densit y maximum a possi ble position to trap 
the protoplanets (jMasset et al. 20061 ) . 

Several mechanisms can lead to the density maximum 
in the gas disk, e.g., the inner disk cavity. At the place 
where the toque induced by the stellar magnetic field 
dominates over disk internal stresses, material in the disk 
will fall o nto the surfa ce of the star thus the disk was 
truncated ()Konigl 199li ) . This place is roughly the coro- 
tation radius of the star, with a maximum stellar distance 
~ 9 stellar radii. Before the protostar moving to the main 
sequence, the radius is 2-3 times larger. Thus the inner 
disk truncation of a main sequence star with mass 0.5 
Mq is about 0.065 AU to 0.085 AU. In our simulations, 
we choose 0.1 AU as the boundary of the inner cavity of 
the gas disk. Once embryos reach this boundary, type I 
migration will be stopped effectively. 

Near the inner edge of the dead zone {ocrit) is an- 
other location that induces density maximum in the 
gas disk. During the classical T Tauri star (CTTS) 
phase, inside flcrit, the protostellar disk is thermally 
ionized, so it is totally active zone. While outside 
O'crit, stellar X-rays and diffuse cosmic rays ionize a 
surface layer with a thickness = lOOg cm~^ on either 
side of the protostellar disk (active layer). Between 
the active layers is the dead zone nea r the midplane 
(jBalbus fc Hawlev 19911 : iGammie 19961 ). Near the in- 
ner edge of the dead zone in the midpla ne, a pos- 
itive density grad ient is induced (Krct ke fc Lin 20071 : 
iKretke et al. 20091 ). For the protostellar d isk in ad 
hoc a-prescription (jShakura fc Sunyaev 19731 ) , the mass 
accretion rate is constant across the disk, Mg — 
ZiraCshTig / , where Cs and h = Cs/^k are the sound speed 
of the midplane and the disk scale-height, respectively. 
As the value of a increases from dead zone (~ 0.006) to 
active zone (~ 0.018) (Sano et al. 2000), the gas density 
in the midplane of the disk increases from active zone 
to dead zone about two times. Then a density maxi- 
mum appears in the inner edge of the dead zone which is 
helpful to halt the embryos undergoing type I migration. 

Due to the density profile enhancement at the inner 
edge of magnetic inactive zone (acrit 0.2 — 0.3 AU), 
embryos undergoing type I migration are expected to be 
stalled there, which makes the location around flcrit an 
ideal place for the accumulation of embryos. Cohesive 
mergers among embryos may finally form terrestrial plan- 
ets. As the location of Ocrit around dwarf stars with mass 
~ 0.5 Mq is roughly in the HZ, the forming planets may 
be most possibly habitable. 

In this paper we study the formation of habitable 
planets due to the above mechanism with a model sys- 
tem OGLE-06-109L. Among the planets detected by mi- 
crolensing, OGLE-06-109L system is the first one with 
observed multiple planets. It is 1490 pc away from the 
Sun, with a star of ^ 0.51 Mq (solar mass) and two 
planets of 0.71 Mj (Jupiter mass) and 0.27 Mj in the 
orbits of 2.3 AU and 4.6 AU, respectively (Gaudi et al. 
2008). The eccentricity of the out planet is 0.11 deduced 
from observation. In the first paper (Wang et al. 2009, 
hereafter WZZ09), we have studied the eccentricity for- 
mation of the system. The planetary scattering model 



or the mean-motion resonance crossing model can be ac- 
count for the eccentricity formation. According to our 
models, the final eccentricity of inner planet (cf,) may 
oscillate between [0-0.06], comparable to that of Jupiter 
(0.03-0.06). 

According to Kasting et al. 1993, HZ is estimated to 
be [0.25-0.36 AU] or larger in OGLE-06-109L system. In 
WZZ09, we have demonstrated orbits in center of HZ 
with semimajor axis a £ [0.28 AU, 0.32 AU] would be 
out of HZ due to the high eccentricities excited by the 
secular resonance of two giant planets in the system. 
An additional inner planet or outer planet would sup- 
press the eccentricity perturbation and greatly improve 
the prospects for habitability of the system (Malhotra & 
Minton 2008). Orbits in the rest region of HZ remain 
in HZ up to 10 Myr integrations. Based on the results, 
OGLE-06-109L is a hopeful candidate system for host- 
ing habitable terrestrial planets. The organizing of the 
paper is as follows. In section 2 we present the model, 
method and initial setup of the simulations. The simu- 
lation results are shown in Section 3. Section 4 presents 
our conclusions. 

2. MODEL 

2.1. Disk Model 

To model the effect of density increment near Ocrit in 
the gas disk, we formulate the gas disk as a viscous one 
with constant mass accreting to the s tar. Its surfa ce 
density at stellar distance a is given as (iPrinde 19811) 



M 



'i'Kv{a) ' 



(1) 



where M is the stellar accretion rate, v{a) is the effective 
viscosity. According to the results from observation data 
of the young cluster p-Oph, the star accretion rate M 
can be estimated to (Natta 2006) 



M^4x lO"8(^)i-8M0yr-i. 

^Mq' °^ 



(2) 



So we adopt a baseline M = 1 x 10"^ Moyr"^ for OGLE- 
06-109L system. However, during the earlier stage of T 
Tauri stage, M c ould be 1-2 order of magnitude larger 
(|Hartmann 20071 ). 

The effective viscosity z/(a) = acsh, where a, c^, and 
h = Cs/^ refer to the efficiency factor of angular momen- 
tum transport, sound speed at the mid plane, and the 
isothermal density scale height, respectively (Shakura & 
Sunyaev 1973). To model this effect, we let aum and 
ttdcad denote the a-values of the MRI active and dead 
regions, respectively. The effective a for the disk is mod- 
eled as (Kretke & Lin 2007; Kretke et al. 2009) 



acff (a) 
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^crit 
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) + l]+aMRi, (3) 



where erf is the error function, O.lCcrit is thought as the 
width of the transition region. In this paper, we adopt 
aMRi = 0.01,adoad = 0.001 (Sano et al. 2000). In ah 
the simulations, we think the embryo to be run into the 
central star if its semimajor axis is less than 0.1 AU. 

The location of the inner edge of the MRI dead zone, 
acrit, varies with the disk temperature, kinematics and 
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Fig. 1. — The density profile of the disk with different star 
accretion rate (M). In the region outside the density maximum, the 
density profile satisfies Eg oc r~^/^. The blue solid line represents 
Af = 1 X 10~* MQyr~^, the location of the density maximum is 
at 0.22 AU. The green dot line shows the density profile of M = 
5 X 10~* MQj/r , the density maximum location is at 0.3 AU. 
The red dash line means M = 2.5 X lO"'^ MQyr~^ , the location of 
the density maximum is at 0.62 AU. The inner edge of the disk is 
in the location of 0.1 AU. 



mass accretion rate, etc. Here we adopt the expression 
from Kretke et al. (2009), 



Ocrit = 0.16 AU( — 



^ N4/9/ ^* Nl/S/ Q'MRK -i/s 



lO-8M0yr- 



0.02 



. (4) 

During the evolution of T-Tauri star and its disk, M will 
decrease on average. So the location of 

^crit 1^ migrating 

inward as the time proceeds. 

Combing with the empirical minimum mass solar neb- 
ula (hereafter MMSN, Hayashi 1981), we obtain the fol- 
lowing surface density for the circumstellar disk (shown 
in Figure 1): 



S.-So(^) (^ 
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where Sq = 883 g cm ^ is the initial density correspond- 
ing to that at 1 Myr for M = 5 x 10"® A/gyr^^ due 
to disk accretion, photoevaporation or planet formation, 
the gas is depleted with a timescale r^ep ~ 3 Myr (Haisch 
et al 2001), t is the evolution time. 

2.2. type I migration and e- damping 

For a planet embedded in a geometrically thin 
and locally isothermal disk, angular momentum ex- 
changes between the planet and the gas disk will 
cause a net momentum loss on the planet, which 
results in a fast and so-called type I i nigration of 
the planet (iGoldreich fc Tremaine 19791 : IWard 19971: 
iTanaka et al. 2003 ) . Some mechanisms are proposed re- 
cently to reduce the speed or even reverse the direction 
of migration. Laughlin et al. (2004), and Nelson & Pa- 
paloizou (2004) proposed that, in the locations where 
the magnetorotational instability (MRI) is active, grav- 
itational torques arising from magnetohydrodynamical 
turbulence will contribute a random walk component to 



the migratory evolution of the planets, thus prolong the 
drift timescale. Analytical estimation through Fokker 
Planck approach indicates the survival rate of planet 
under type I and stochastic torques is around 1 — 10% 
(Adams & Bloch 2009). Paardekooper & Mellema (2006) 
noticed that the inclusion of radiative transfer can cause 
a strong reduction in the migration speed. Subsequent 
investigations (Baruteau & Masset 2008; Kley & Crida 
2008; Paardekooper & Papaloizou 2008) indicated that 
the migration process can be slowed down or even re- 
versed for sufficiently low mass planets. Through full 
three-dimensional hydrodynamical simulations of embed- 
ded planets in viscous, radiative disks, Kley et al. (2009) 
confirmed that the migration can be directed outwards 
up to planet masses of about 33 . 

Taken these uncertainties into considerations, we adopt 
a standard description of type I migration with a speed 
reduction factor fi . The timescale of t ype I migration 
for an embryo m is (|Tanaka et al. 2002D 
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where r, e, h, f2, and Sg are the stellar distance, eccen- 
tricity of the embryo, scale height of the disk, the Kepler 
angular velocity, and the surface density profile of gas in 
the disk, respectively, /3 — — dlnSg/dlna. According to 
the density profile in equations ([3]) and ([5]), the factor 
/3 will be negative when crossing the location of density 
maximum to the inner disk. So when the embryo goes 
across the density maximum, the embryo may be sus- 
tained in the maximum location. In our simulations, we 
consider /i is in the range of [0.001,1]. The type I migra- 
tion effect is incorporated with an additional acceleration 
to any embryo with mass < 10 (Cresswell & Nelson 
2006), 

V 

Fmig/ = , (7) 

where v is the velocity vector of m in the stellar-centric 
coordinated] 

Interactions between embryos and the gas disk may 
damp the eccentricities of the embryos (Goldreich & 
Tremaine 1980). The timescale of the eccentricity damp- 
ing for an embryo with mass m < 10 can be de- 
scribed as (Cresswell & Nelson 2006), 



( . )cdap 
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where r, e, h, are the distance from the star, eccen- 
tricity of the embryo, scale height of the disk and the 
Kepler angular velocity, respectively, Qe = 0.1 is a nor- 
malization factor to fit with hydrodynamical simulations. 
This effect is included by add an extra-force term to the 
equations of embory motion: 
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2.3. 



Initial configurations of the system 

We assume two giant planets (m;,, to^) have already 
formed with observed masses of OGLE-09-106L b, c ini- 
tially at 4 AU and 8 AU, respectively. The giant planets 
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undergo type II migration with a fixed timescale about 1 
Myr (WZZ09). An additional embryo with mass 16.8 M© 
is put between them. A planetary scatter will occur 
between one giant planet and the additional embryo at 
about 5 Myr, just as in WZZ09, to account for the ob- 
served eccentricities of the two giant planets. 

To model the formation of planets in the terrestrial 
region, we put 18 embryos with isolation masses (Ida & 
Lin 2004) 



Mi, 



(^^)3/4(^^-l/2^^ (10) 



in the inner region, where 7ico is the volatile enhancement 
with a value of 4.2 or 1 for material exterior or interior 
to the snow line (0.68 AU for OGLE-06-109L system), 
respectively. The embryos with the isolation mass are 
located in 10 mutual Hill radii from each other. For a 
dwarf star with mass 0.51 Mq, we assume a disk with 
fd < 2 . For fd — 1, the masses of embryos range from 
0.084 Me to 3.33 M® at [0.39 AU, 3 AU], with a total 
mass of 19.1 M^. For fd = 2, the masses of embryos 
range from 0.175 M® to 9.42 M® at [0.26 AU, 3 AU], 
with a total mass of 38.5 Mq. The outmost embryo 
in the inner region is in the orbit 3.5 Hill Radii away 
from mf,. The giant planets and embryos are all in near- 
coplanar and near-circular orbits initially (with initial 
eccentricities e — 10~^ and inclinations i — e/2). 

2.4. Perturbation from two giant planets 

The effect of giant planets is twofold: for embryos in 
nearby orbits, perturbations from two giant planets will 
excite their eccentricities, which may be helpful for the 
merge of small embryos into ones. The maximum eccen- 
tricity that an embryo initially in circular orbits exited 
under the perturbation of out planet is (Mardling 2007) 

{5/2){ae/ac)ec€-^ ^^^^ 



1 - y^ae/ac{me/mc)e 

where subscribes of e and c denote those of the em- 
bryo in inner orbit and the giant planet, respectively, 
e = a/1 — e^. According to equation dTTI) . the excitation 
of initial embryos located from [0.26, 3] AU is [0.01, 0.10] 
from TOb and [0.007, 0.08] from nic- Thus the effect of 
inner giant planet is non-negligible for embryos around 3 
AU. 

Another effect that may contribute to the formation of 
planets in habitable region is the secular resonance. Due 
to the presence of gas disk and mutual perturbation, the 
orbits of two giant planets will precess with eigenfrequen- 
cies (si,S2) (Heppenheimer 1980), 

A22)± [(An -^22)^ +4^12^21 



Sl,2 



where 
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All = 2nbabmcNbc - (7rF/4ab)cro, 
A12 = -2nbabmcPbc, 
A22 = 2ncacmbNbc - [ttF / Aac)(Jo, 
A21 = -2ncacmbPbc, 

where cto — So(aeff/10~'^)~^ exp(— t/rdep). We define 
a_ — min(ab,ac) and a+ = max(af,,ac), so that 

Nbc = ^{a^/al)b''^)\{a^/a+), 

Pbc = lia-/a\)bfl2{a^/a+), 
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Fig. 2. — Location of secular resonance in the inner region of 
the system at different time during the evolution, (a) The precess 
eigenfrequencies si, S2 of the two giant planets due to different 
density of gas disk So and mutual perturbation, (b) At t=0, the 
proper frequency of test particle g compared to si, S2- There is 
no point of intersection between g and si, S2. Secular resonance 
cannot happen in the inner region of the system, (c) At t=3.76 
Myr, the proper frequency of test particle g compared to si, S2- g 
is close to si. (d) At t=4.48 Myr and 10 Myr, the proper frequency 
of test particle g compared to si, S2- Secular resonance happens 
at the location of ~ 0.25 AU. 



where b\^j^ is the Laplace coefficients. F = 4.377345 for 

surface density slope /3 — 3/2. 

The proper frequency of test particle associated with 
the longitude of perihelion is 



2npap ■ {nibNpb + nicNpc) 



{TTF/4ap)ao, (13) 



where subscript p refers to the test particle, terms with 
Npb, Npc are due to the perturbation of two planets, term 
with (To is due to the gas disk. The definition of Npb and 
Npc are similar to Nbc- When the proper frequency of 
an embryo meets with either of the two eigenfrequencies, 
i.e., g = Si, (i=l, 2), secular resonance occurs, which 
may excite the eccentricities of embryos. 

The two eigenfrequencies (si, S2) are plotted in Figure 
2a. As we can see, the orbits of two planets precess ret- 
rogradely (s^ < 0) unless the gas disk is almost depleted 
(do < 10 gcm~^), so the procession rate is slow (with 
period > 2000 yr) and prograde. Also the precession fre- 
quency g located inside 2.0 AU dose not equal to si, S2 
unless Sg < 20 gcm~^ (Figure 2b-d). This occurs when 
i > 3.8 Myr. At i = 10 Myr, the gas disk is effectively 
depleted, either si or S2 may meet g at some locations of 
habitable zone. As the secular resonances at habitable 
zone occurred only at the epoch when the gas disk is al- 
most depleted, we omit the effect of self gravity of gas 
disk in the following simulations. 
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3. SIMULATIONS AND RESULTS 

We numerically integrate the full equations of planet 
motions with a time-symmetric Hermit scheme (Aarseth 
2003). Regularization technique is used to handle the 
collision between embryos, so all the embryos have their 
physical radii. As long as the mutual distance between 
two embryos is smaller than the sum of their physical 
radius, we assume they are merged. 

We perform four Groups of simulations with differ- 
ent star accretion rate (M) and the heavy element scal- 
ing factor with respect to the MMSN model (/d), see 
Equation pO| . In our simulations we choose M = 
1 X 10"*^, 5 X 10~®, and 2.5 x 10"'^ Mq/jv, correspond- 
ing to different stages of T Tauri stars so that higher M 
corresponds to earlier stage. The conditions in different 
Groups are listed as follows: 

Group 1 : /<j 2; M = 5 x IQ-^Me/yr responds to the 
green dot line in figure 1. 

Group 2 : /d = 2;M = lx IQ-^Mg/yr responds to the 
blue solid line in figure I. 

Group 3 : /d = 2; M = 2.5 x lO^'^Me/yr responds to 
the red dash line in figure 1. 

Group 4:/d = l;M = 5x lO^^Mg/yr responds to the 
green dot line in figure 1. 

We perform 31 runs with the same disk model except 
varying /i from 0.001 to 1 in each Group. The results 
are presented as follows. 

3.1. Formation and orbital configurations of terrestrial 

planets 

One of the most important factors that affects the for- 
mation of planets in habitable zone is the speed of type 
I migration reduced by /i. If the speed is very slow, 
the embryos will not migrate too much. Fig. 3 shows a 
typical run of Group 1 with very low migration speed 
(/i = 0.004). Planetary mergers mainly occur due to 
scattering effect of embryo inside two giant plants at 
t ^ b Myr (WZZ09). As we know from the density pro- 
file in Figure 1, when M = 5 x 10~* M0/yr, the location 
of density maximum is at 0.3 AU. At this density max- 
imum, the type I migration of embryos can be stalled 
due to the negative /? in equations ([S]). The inward mi- 
gration of a second planet may trap the stalled planet 
into its mean motion resonance, which might increase 
the eccentricities of both planets. As we can see, most 
of the planets survived at i = 10 Myr with a < 1 AU are 
trapped into first order resonance (mainly 4:3). Among 
them, one is in the HZ and two of them are at the edge 
of the HZ in this case. 

With the increase of the type I migration speed, less 
and less embryos are survived or trapped in HZ before 
the gas disk was almost depleted. For example, Figure 4 
shows a typical run with migration reduced factor /i = 
0.08. Only two terrestrial planets survived at the inner 
region (a < 1 AU). Interestingly, they are in 2:1 MMR. 
The inner planet with mass (7.5 M^) is in the HZ. As 
the migration speed increases again, /i — 0.25, only one 
planet survived in the HZ (Figure 5a). It is stalled under 



type I migration at the inner edge of MRI dead zone 
(Figure 5b). 

According to Figure 3b, many planet pairs may be 
trapped in MMRs under migration. Table 1 shows the 
statistics of the final configurations for the simulations in 
the four Groups. We classify the system being in either 
MMRs, aligned apsidal resonance (AARs) or anti-aligned 
apsidal resonance (anti-AARs) if there is at least one cou- 
ple of planets trapped into MMRs, AARs or anti-AARs. 
Table 1 records the probability of a planet system that 
is in MMRs, AARs or anti-AARs in the 31 runs of each 
group (24 runs in Group 3). As we can see, more than a 
quarter of the systems contain at least a planet pair that 
is in MMRs, AARs or anti-AARs. 

3.2. Stable time of the system 

For systems with only one terrestrial planet, we can 
get the stable time from the results in WZZ09. Because 
the distance of the single terrestrial planet to the cen- 
tral star is closer than 0.56 AU in our results. All the 
single terrestrial planets systems are stable in more than 
10* years. Based on Zhou et al. (2007) about the cross- 
ing time of planet system, we estimate the stable time 
of multiple terrestrial planet system at the end of the 
simulations. We test each couple of the planets in the 
systems and use the average eccentricity, average mass 
and relative Hill radius for the two planets next to each 
other in the estimation. Results are showed in table 2. 
More than 60% couples of the planets in the system will 
survive in more than 10* years. According to the results 
of the stable time of the planets, we get the stable time 
of the planetary system. We define the shortest time of 
the planet couples surviving in the system as the stable 
time of the system. The results are also showed in the 
last column of table 2. More than 41% system will be 
stable in more than 10* years. 

3.3. Water contents of terrestrial planets 

For terrestrial planets in habitable zone, one important 
parameter is the water content inside the rock. To model 
the distribution of water contents within embryos, we ini- 
tially assume the inner disk is water-poor, embryo within 
0.5 AU is 0.001% by mass, while embryo beyond 0.625 
AU is wet with the water content 5%, embryo between 
them is moderately wet with 0.1% water. We further 
assume, the water content will be decreased by 30% dur- 
ing each collision (Melosh 2003). Thus the migration and 
merge of embryos may re-distribute the water contents. 

Figure 6 shows the simulation results of Group 1. The 
color map represents the water distribution by mass 
(Raymond et al. 2004). Panel (al) and (bl) show the 
distribution of the habitable planets changed with /i . If 
type I migration of the embryos are fast, /i > 0.2, most 
of the small embryos originally in the HZ were run into 
the central star, the only survival planet with high mass 
(^ 7 Mq) was migrated from outside (> 2.2 AU). So 
they have relative high water contents. When /i < 0.07, 
the slow migration of embryos are easily to be perturbed 
by either the secular perturbations or secular resonance 
induced by the two giants. So planetary merges occurred 
occasionally, and the motion of the embryos are chaotic. 
As planetary merges reduce the water contents of the 
rocks, such formed planets with masses > 7M(q are drier 
than neighbor embryos . 
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Fig. 3. — A typical run with /i = 0.004 in Group 1. There are nine terrestrial planets with masses less than 7M^ left in the system. 
One is in the HZ and two of them are at the edge of the HZ. At the end of the simulation they trapped into 4:3 and 3:2 MMRs respectively, 
(a) Evolution of the semi-major axes. The gray band represents the HZ and the two blue dashed lines indicate the location of secular 
resonance caused by the giant planets, (b) The final locations of the terrestrial planets in the system. The ordinate of each planet is chosen 
to fit the gas density profile for illustrative purposes. MMRs are labeled. 




Fig. 4. — A typical run with /i = 0.08 in Groupl. There are two terrestrial planets left in the system. One is in the HZ. At the end 
of the simulation they trapped into 2:1. (a) Evolution of the semi-major axes. The gray band represents the HZ and the blue dashed line 
indicates the location of secular resonance caused by the two giant planets and the other terrestrial planet, (b) The final locations of the 
terrestrial planets in the system. The ordinate of each planet is chosen to fit the gas density profile for illustrative purposes. 



TABLE 1 

Statistic of the final configurations of the simulations in the four groups. f-Pb and f-Pc show the final average 
locations and eccentricities of giant planet b and c respectively. f-mmrs is the probabilities of planets (in all planets) 

AND THE SYSTEM (wiTH AT LEAST ONE PLANET PAIR) CAPTURED INTO MMRS. We CONSIDER THE SYSTEM IS IN MMRS IF THERE IS AT 
LEAST ONE COUPLE OF PLANETS TRAPPED INTO MMRS. F-AAR AND F-ANTlAAR ARE THE STATISTIC OF THE PLANETS AND SYSTEMS 
TRAPPED INTO ALIGNED APSIDAL RESONANCE AND ANTIALIGNED APSIDAL RESONANCE IN THE FINAL CONFIGURATIONS. 



ID 


f-Pb 


f-Pc 


f-MMRs 


f-AAR 


f-antiAAR 




a (AU), e 


a (AU), e 


Planet, System 


Planet, System 


Planet, System 


Gl 


2.474, 0.076 


5.267, 0.108 


63.4%, 17/31 


62.7%, 18/31 


56.2%, 15/31 


G2 


2.547, 0.087 


5.272, 0.132 


47.6%, 16/31 


48.4%, 12/31 


48.4%, 16/31 


G3 


2.486, 0.065 


5.166, 0.101 


30.8%, 7/24 


46.2%, 6/24 


27.4%, 6/24 


G4 


2.000, 0.087 


5.040, 0.107 


23.0%, 11/31 


56.0%, 13/31 


41.0%, 13/31 
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Fig. 5. — A typical run with /i = 0.25 in Groupl. There is one terrestrial planet left in the system and it is in the habitable zone, 
(a) Evolution of the semi-major axes. The gray band represents the HZ and the two blue dashed lines indication the location of secular 
resonance caused by the two giant planets, (b) The final location of the terrestrial planet in the system. The ordinate of the planet is 
chosen to fit the gas density profile for illustrative purposes. 



TABLE 2 

Statistic of the stable time of the multiple terrestrial planet systems. The first three columns means the statistic of 

THE planet couples AND THE LAST COLUMN REPRESENTS THE RESULTS OF THE PLANETARY SYSTEMS IN EACH GROUP. 



ID 


> 10" yr 


> 10' yr 


> 10« yr 


system stable time > 10** yr 


Gl 


92.6% 


84.4% 


63.9% 


41.67% 


G2 


90.4% 


86.2% 


69.1% 


47.82% 


G3 


91.1% 


83.3% 


77.8% 


50% 


G4 


92.8% 


92.8% 


87% 


75.86% 



water content % 



0.40 
r 0.36 
S 0.32 
0.28 
0.24 



S 

£ 8 
S 

4 



(a1) 







I I 1 1 1 1 1 1| 



I I 1 1 1 1 III 



I I I 1 1 1 III 



1E-3 



0.01 



0.1 





Fig. 6. — The results of Group 1 with M = 5 X 10~* MQyr~^ at the end of our simulations (t=10 Myr). (al): The semimajor axes 
of the habitable planets distribution changed with /i. (bl): The Mass of the habitable planets distribution. (a2): The number of the 
terrestrial planets survived in the system. (b2) The average semimajor axes of terrestrial planets distribution changed with /i . (c2) The 
average masses of the terrestrial planets distribution. The color map in the left panels represents the water distribution by mass. The error 
bars in Panel (al) and (b2) show a(l-e) and a{l+e). 
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Panels (a2), (b2), and (c2) of Figure 6 are the dis- 
tribution of the terrestrial planets. From Figure 6a2, if 
the type I migration speed is slow, there will be many 
embryos left in the inner region, fi ^ 0.1 is the crit- 
ical value fieri above which only one or two terrestrial 
planets formed and survived in the system. The average 
semi-major axes of the terrestrial planets showed in Panel 
(b2) are small when /i > 0.1, within the snowline (0.68 
AU in OGLE-06-109L system), so most of the terrestrial 
planets keep in the inner region of the system. As it's 
possible to form planets with mass > 10 M^, they might 
have the chance to accretion gas and form Neptune-size 
planets if there is still enough gas in the disk. 

In Group 2 and 3, we change the star accretion rate to 
M = 1 X 10"* Mg/yr and M = 2.5 x 10""^ Mg/yr respec- 
tively. This corresponds to the formation of terrestrial 
planet occurred either later or earlier so that the mass ac- 
cretion rate of the star is either smaller or larger. Figure 
7 shows the results of the two groups. When /i > 0.1 of 
Group 2, there are one or two planets left in the system, 
and the inner most terrestrial planet with mass ^ 9.4 Af^ 
locates around 0.22 AU, the density profile maximum 
around the HZ. For Group 3, M = 2.5 x lO-'^Mg/yr, 
the location of density maximum is 0.62 AU. There is no 
terrestrial planet left in the system for /i > 0.2. 

From Figure 7c 1 and c2, we can get that if /i < 0.1, 
there are more than one terrestrial planets survival in the 
system. It also agrees with the result of Terquem & Pa- 
paloizou (2007) that the hot super-Earths or Neptunes 
formed by mergers of inwardly migrating cores are most 
likely not isolated. The average semimajor axes in Fig- 
ure 7dl and d2 show that most of the terrestrial planets 
are distributed inside the snowline (0.68 AU). The water 
content in the region of /i < 0.1 is lower than that of 
/i > 0.1, in which case the terrestrial planets are mainly 
from outside orbits with high water contents. 

Figure 8 shows the results of Group 4 for the embryos 
with /d = 1. Because fd is related to the mass of the 
embryos initially, the mass of the habitable planets will 
be smaller than other Groups. From the timescale of 
type I migration, equation (6), the migration speed is in 
proportion to /i x m. So when /d = 1, the migration 
is slower than the runs in other Groups. In Group 4, if 
/i is less than 0.2, most of the habitable planets left in 
the system are merged from small embryos for several 
times. The water contents of the planets are low. When 
fd ~ 1, the initial masses of all embryos are less than 
3.3 Mq. At last the average mass of Group 4 is about 
4 Mq. The habitable planet whose mass higher than 7 
is formed from merge with low water content. In 
this Group, there are more than one terrestrial planets 
survival in the range of /i e [0.001, 1]. 

4. CONCLUSIONS AND DISCUSSIONS 

In this paper, we investigate the formation of terrestrial 
planets under the perturbation of giant planets around 
dwarf stars, with a model system OGLE-06-109L. We as- 
sume two giant planets already formed in outside orbits. 
The embryos are initially embedded in an exponentially 
decaying gas disk and will under type I migration with 
a speed reducing factor fi over the standard (linear) es- 
timation. The presence of giant planets may excite the 
eccentricities of embryos in nearby orbits through secular 
perturbations, and induce secular resonance in the hab- 



itable zone regions when the gas disk is almost depleted. 

The speed of type I migration plays an important 
role in determine the number and orbital architecture 
of terrestrial planets. When /i < 0.1 (0.1 times of the 
standard speed), there may be several terrestrial planets 
formed in the inner region (< 1 AU). Massive terrestrial 
planets formed by mergers of smaller mass embryos, thus 
they have lower water contents compared with neighbor- 
ing embryos with small masses. Some of the planet pairs 
are trapped either in MMRs, or AARs or anti-AARs. 
The formation of planets in MMR is due to the halt of 
inner planets at the inner edge of MRI dead zone, or due 
to the convergent migration during the evolution. Ac- 
cording to our simulations, most of the formed MMRs 
are first order resonances which have wider resonance re- 
gions. In the case that type I migration speed is fast 
(/i > 0.1), there will be very few (or sometimes no) 
terrestrial planets in the inner region. They are mainly 
migrated from outer regions so they have higher water 
contents as compared with neighboring embryos. Thus 
through the water contents of the terrestrial planets, we 
can deduce their origin. 

The change of star accretion rate refiects the location 
of the density maximum, which affects the survival of 
the planets in the habitable region [0.25, 0.36] AU in 
OGLE-06-109L system. If the giant embryos in the inner 
region (< lAU) of the system formed in the earlier (for 
example. Figure 7a2-e2, M = 2.5 x 10"'^ MQ/yr), the 
surface density of the disk is still high so that the high 
speed type I migration (with /i > 0.1) may drive all the 
embryos into the star. Only lower speed migration and 
later formed embryos are easier to retained. However, 
the model of type I migration described by equations 
(O and (O are linear. When more realized model is 
considered, e.g., considering the outward migration in 
viscous, radiative disks (Kley et al. 2009), things will be 
more complicated and interesting. 

Another important factor is the initial conditions of 
the embryos. In our simulations, the embryos we used 
are in isolation masses. In order to test the reasonability 
of the initial conditions, we run with more embryos with 
smaller initial masses using /i = 0.004. Figure 9 shows 
the result. Finally, we get a planet with mass of 4.4 M^ 
in the habitable zone which is comparable to the result in 
figure 3. And the system will be stable in more than 10* 
years. We test other runs and find that the results are 
similar to the run when /i = 0.004. So our results with 
the embryos in isolation mass initially are reasonable. 

According to our simulations, by checking the water 
contents of terrestrial planets in OGLE-06-109L (and 
other similar systems), one can understand their forma- 
tion history, and also restrict the parameters such as /i 
used in the planetary formation theory. For the OGLE- 
06-109L system, the snowline is about 0.68 AU, very close 
to the central star. If one or no terrestrial planet is ob- 
served, basically they might formed through inward mi- 
gration from outside the snow line, thus they will have 
relatively high water contents (^ 0.5%), also the type I 
migration speed is fast (/i > 0.1). If many planets are 
observed in the inner region (< 1 AU), then most prob- 
ably, they are formed through merges of isolation cores, 
thus they will have relatively low water contents ~ 2%. 
The presence of multiple terrestrial planets also indicates 
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Fig. 7. — The results of Group 2 (left panels ) and Group 3 (right panels) at the end of our simulations (t=10 Myr). Panel (al) and 
(bl) are the semimajor axes and mass distribution of the habitable planets in Group 2 with M = 1 X 10~* MQyr~^ . Panel (cl), (dl), and 
(el) are the number, average semimajor axes, and average mass distribution of the terrestrial planets in Group 2. Panel (a2) and (b2) are 
the results of habitable planets in Group 3 with M = 2.5 X 10~^ MQyr~^. Panel (c2), (d2), and (e2) are the results of terrestrial planets 
of Group 3. The error bars in panel (al), (a2), (dl), and (d2) show a(l-e) and a(l+e). 
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Fig. 8. — The results of Group 4 at the end of our simulations (t=10 Myr). Panel (al) and (bl) are the habitable planets distribution 
of Group 4 with fd = 1. Panel (a2), (b2), and (c2) are the results of terrestrial planets in Group 4. The error bars in Panel (al) and (b2) 
show a(l-e) and a(l+e). 



a low speed migration of the embryos (/i < 0.1). 
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